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Abstract 

The Lighthill acoustic analogy, as embodied in the 
Ffowcs Williams-Hawkings (FW-H) equation, is com- 
pared with the Kirchhoff formulation for moving sur- 
faces. A comparison of the two governing equations 
reveals that the main Kirchhoff advantage (namely 
nonlinear flow effects are included in the surface in- 
tegration) is also available to the FW-H method if the 
integration surface used in the FW-H equation is not 
assumed impenetrable. The FW-H equation is ana- 
lytically superior for aeroacoustics because it is based 
upon the conservation laws of fluid mechanics rather 
than the wave equation. This means that the FW-H 
equation is valid even if the integration surface is in 
the nonlinear region. This is demonstrated numeri- 
cally in the paper. The Kirchhoff approach can lead 
to substantial errors if the integration surface is not 
positioned in the linear region. These errors may be 
hard to identify. Finally, new metrics based on the 
Sobolev norm are introduced which may be used to 
compare input data for both quadrupole noise calcu- 
lations and Kirchhoff noise predictions. 

Introduction 

A great deal of progress has been made in recent 
years toward the prediction of rotating-blade noise 
through methods utilizing first principles. Several rea- 
sons account for this progress. First, a detailed and 
fundamental understanding of how rotor blades gen- 
erate noise has been gained through several acoustic 
wind-tunnel and flight tests. Secondly, a rigorous the- 
oretical basis for predicting noise generated by rotating 
blades has been developed. In fact, several prediction 
methodologies with a solid physical and mathemati- 
cal basis are currently available: formulations based 
upon the Lighthill acoustic analogy 1 (in particular the 
Ffowcs Williams and Hawkings (FW-H) equation 2 ) 
and Kirchhoff formulations for both subsonic and su- 
personic moving surfaces . 3,4 
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The FW-H equation is an exact rearrangement of 
the continuity equation and the Navier-Stokes equa- 
tions into the form of an inhomogenous wave equa- 
tion with two surface source terms, known as thick- 
ness and loading sources, and a volume source term 
(the quadrupole source from the original Lighthill the- 
ory). Although the quadrupole source contribution is 
insignificant in many subsonic applications, substan- 
tially more computational resources are needed for vol- 
ume integration when the quadrupole source required. 
The Kirchhoff formulation for moving bodies is also an 
inhomogenous wave equation with source terms dis- 
tributed on a surface which encloses all of the physical 
sources. The Kirchhoff formulation is attractive be- 
cause no volume integration is necessary. 

Although it is useful to have more than one formula- 
tion available to predict noise, there is no clear consen- 
sus of which to choose for a particular application. A 
recent numerical comparison by Brentner et al . 5 of the 
helicopter rotor noise prediction code WOPWOP +, 6-8 
which uses a FW-H based formulation including an 
approximate quadrupole calculation, with a rotating 
Kirchhoff code RKIR 9,10 has shown that both meth- 
ods can predict the rotor noise equally well. In that 
work, however, neither method was demonstrated to 
be clearly superior. 

The main purpose of this paper is to analytically 
compare these two acoustic prediction methodologies 
and reduce the confusion that currently exists about 
the relationship between the two methods. This in- 
cludes a comparison of how the governing equations 
are derived, highlighting the differences in the deriva- 
tions. Both analytical comparison and numerical 
comparisons are necessary to determine whether one 
method has an advantage in terms of efficiency, ac- 
curacy, and robustness over the other. An alternate 
implementation of the FW-H equation is presented 
which combines the advantages of both the traditional 
formulation of the FW-H equation and the Kirchhoff 
formulation. Finally, a useful metric for comparing 
formulations will be outlined. 
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Advantages and Disadvantages 

It is important to briefly consider the advantages 
and disadvantages of both the FW-H and Kirchhoff 
formulations at this point in order to understand the 
motivation for a more in depth analysis. 

FW— H Description 

The FW-H approach has several advantages over 
the Kirchhoff method. First, the three source terms in 
the FW-H equation each have physical meaning which 
is helpful in understanding the noise generation. The 
thickness noise is determined completely by the geom- 
etry and kinematics of the body. The loading noise is 
generated by the force acting on the fluid due to the 
presence of the body. The classification of thickness 
and loading noise is related to the thickness and load- 
ing problems of linearized aerodynamics. Thus, this 
terminology is consistent with that of aerodynamics. 
The quadrupole source term accounts for nonlinear ef- 
fects (i.e., nonlinear wave propagation and steepening, 
variations in the local sound speed, noise generated 
by shocks, vorticity, and turbulence in the flow field, 
etc.) 11-13 

All three source terms are interdependent, yet their 
physical basis provides information to design quieter 
rotors. The separation of source terms also is an ad- 
vantage numerically because not all terms must be 
computed at all times if it is known that a particu- 
lar source does not contribute to the sound field (e.g., 
for low-speed flow the quadrupole may be neglected, 
in the rotor plane thickness noise is dominant, etc.). 
A final advantage of FW-H based formulations is that 
these formulations are relatively mature and have ro- 
bust numerical algorithms. The main disadvantage of 
the FW-H method is that to predict the noise of bod- 
ies moving at transonic speeds the quadrupole source 
must be included. This is a disadvantage because the 
quadrupole — which is a volume source — ultimately re- 
quires a volume integration of the entire source region. 
Volume integration is computationally expensive and 
can be difficult to implement. Although the computa- 
tional effort can be reduced by approximation of the 
quadrupole, 7,8 it cannot be avoided completely. 

Kirchhoff Description 

The Kirchhoff approach does not suffer from this pit- 
fall because it only has surface source terms. Hence, 
the Kirchhoff method has been used for the past sev- 
eral years for the prediction of transonic rotor noise. 
Unlike the FW-H source terms, however, the Kirch- 
hoff source terms are not easily related to thickness, 
loading, nonlinear effects, or indeed any physical mech- 
anisms. They provide little guidance for design. An- 
other disadvantage of the Kirchhoff method is that the 
source surface (Kirchhoff surface) must be chosen to be 


in the linear flow region , such that the input acoustic 
pressure p 1 = p — p 0 and its derivatives dp 1 /dt and 
dp 1 /dn are compatible with linear wave propagation. 
The location of the linear region is not well defined and 
is problem dependent. It would be desirable to place 
the Kirchhoff surface well away from the source region, 
but typically CFD solutions are not as well resolved or 
as accurate away from the body. Hence, the placement 
of the Kirchhoff surface is usually a compromise. 


Analytical Comparison 

Now that the general characteristics of both the 
FW-H and Kirchhoff formulations have been de- 
scribed, a more detailed comparison will be helpful. 
First, we shall consider the development of the gov- 
erning equations of both approaches to gain insight 
into the validity of each type of formulation. Then 
an assessment of an integral formulation for subsonic 
source motion will be considered. 


Governing Equations 


FW-H Equation 

The FW-H equation 2 is the most general form of the 
Lighthill acoustic analogy and is appropriate for pre- 
dicting the noise generated by the complex motion of 
helicopter rotors. The FW-H equation may be derived 
by embedding the exterior flow problem into a prob- 
lem in unbounded space by using generalized functions 
to describe the flow field. To do this, consider a mov- 
ing surface /(x, t) = 0 with a stationary fluid outside. 
The surface f — 0 is defined such that V/ = n, where 
n is a unit normal vector pointing into the fluid. In- 
side f — 0 the generalized flow variables are defined to 
have their freestream values, i.e., 
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where the tilde indicates that the variable is a gener- 
alized function defined throughout all space. On the 
right hand side p, pi^, and Pij are the density, momen- 
tum, and compressive stress tensor, respectively. Note 
that we have absorbed the constant —p 0 &ij into the 
definition of P \j for convenience, hence, for an inviscid 
fluid, P^ = p'Sij. Freestream quantities are indicated 
by the subscript o and Sij is the Kronecker delta. 

Using definitions (1)— (3) , a generalized continuity 
equation can be written 
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where the bar over the derivative operators indicate 
that generalized differentiation (i.e. , differentiation of 
generalized functions) is implied and p' = p — p Q . Also 
note that df/dt = — i? n , df/dxi = hi and 5(f) is 
the Dirac delta function. This generalized continuity 
equation is valid for the entire space — both inside and 
outside of the body. The generalized momentum equa- 
tion can be written 


FW-H equation. The difference is that the domain 
is now considered in terms of wave propagation. The 
surface / = 0 is defined such that all of the acoustic 
sources are contained inside the surface. Then, the 
acoustic pressure p'(x, t) is extended such that 
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Now by taking the time derivative of equation (4) and 
subtracting the divergence of equation (5), followed 
with some rearranging, the FW-H equation may be 
written as the following inhomogeneous wave equation: 

nV<x - <) = 

d 

Ox "F pu{(u n ^n))5(f)\ 

+ 7^[(PoV n + p(u n - V n ))8(f)\ (6) 

where Tij is the Lighthill stess tensor, u n is the fluid 
velocity in the direction normal to the surface / = 0 
and v n is the surface velocity in the direction normal to 
the surface. On the left hand side we use the custom- 
ary notation p' = c 2 p' because the observer location is 
outside of the source region. 

Usually in the derivation of the FW-H equation the 
surface / = 0 is assumed to be coincident with the 
physical body surface and impenetrable (u n = v n ). 
That assumption is not necessary and has not been 
made in equation (6) so that it may be compared more 
directly with the governing equation of the Kirchhoff 
formula for moving surfaces. Ffowcs Williams and 
Hawkings used slightly different mathematical manip- 
ulations, but it is clear from their paper 2 that they 
understood it is not essential to choose the integra- 
tion surface coincidental with the physical body. Re- 
cently di Francesantonio 14 and Pilon and Lyrintzis 15 
have also treated the FW-H on a permeable surface, 
but have used different names to identify the form 
of the FW-H equation given in equation (6). (Pilon 
and Lyrintzis results appear to be incorrect because 
they have substitued p 1 for c 2 p' in some of their source 
terms.) 

Kirchhoff Equation 

The development of the Kirchhoff formulation, due 
to Farassat and Myers, 3 utilizes the same mathemat- 
ical style and rigor as used in the derivation of the 


and the generalized wave equation — which is the 
governing equation for the Kirchhoff formulation — 
becomes 
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where M n = v n /c. In this equation is p 1 must be com- 
patible with the wave equation, hence, equation (8) is 
valid only in the region of the fluid in which the wave 
equation is the appropriate governing equation. 

Source Term Comparison 

It is well known that the wave equation can be de- 
rived directly from the conservation laws of fluid me- 
chanics, but it is our objective in this paper to show 
how equation (8) is related to the FW-H equation, 
equation (6). To that end, we add and subtract terms 
to the inviscid form of equation (6) to manipulate the 
source terms into the form of equation (8). This yields 


□V(x,£) = Q kir + 


B 2 


Oxidx 




o r \ 

- [ PUi(u n - v n )5(f)] + — [pu n 5(f)] . (9) 


If we note that 


d 2 H(f) 

Otdxi 




j^-M(/)) do) 


and utilize the continuity and momentum equations 
we can rewrite equation (9) as 

□V(X,*) = Qkir + ^U-[Tyif(/)] 

+ |b/ - cV]^(/) + | [U - c 2 P ')^S(f)\ 
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This form of the FW-H equation is helpful because 
the source terms that are not found in the Kirchhoff 
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governing equation are easily identified. This is an 
important result of this paper. All of the additional 
source terms are second order and may be neglected 
in the linear flow region. This was precisely Lighthill’s 
original premise — the wave equation is the appropriate 
governing equation outside of a limited source region. 
In fact, when p 1 = c 2 p' equation (11) becomes 

□ VOM) = Qkir + ■ (12) 


Notice that the Heaviside function has been taken out 
of the equation (11) quadrupole source term in the ma- 
nipulations leading to equation (12). The only source 
term remaining which is not in equation (8) is clearly 
second order in the perturbation quantity This 
term would be neglected in the derivation of the wave 
equation from the fluid conservation laws. Hence, we 
have shown that the FW-H and Kirchhoff formula- 
tions are indeed equivalent when the integration sur- 
face for both is placed in the linear region of the flow 
(i.e., where the input data is compatible with the wave 
equation). 

The FW-H equation and the Kirchhoff are quite dif- 
ferent, however, when the integration surface is in the 
source region. The implications of this difference is 
demonstrated later with numerical examples. If the 
FW-H equation integration surface is on the body 
or in the source region, the quadrupole — a volume 
source term — must be included to accurately predict 
the noise. Therefore, we can infer that as we move the 
integration surface of the FW-H equation away from 
the body, the contribution of the volume quadrupole 
contained within the surface must now be accounted 
for by the surface source terms. We shall numerically 
demonstrate this later. 

For completeness, equation (11) can be simplified by 
canceling terms and rearranging. The result is 
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Notice that the surface source terms in equation (13) 
are closely related to equation (8). In fact by sub- 
stituting c 2 p' for p' in the time derivative terms in 
equation (8) and puiUj +p'Sij in the spatial derivative 
terms we can get the surface source terms in equa- 
tion (13). (The momentum equation was used to ex- 
change d(pUiUj + p f Sij)/dxj with —dpui/dt in equa- 
tion (13).) While the correspondence between equa- 
tion (13) and equation (8) is interesting, equation (13) 


has two pitfalls: it is not easily recognized as the 
FW-H equation, and there are no clear connections 
between the form of the source terms and the problem 
physics. 

An Integral Formulation 

Now that the relationship between the FW-H equa- 
tion and the Kirchhoff formulation has been developed 
on the governing equation level, we would like to de- 
velop an applicable integral form which is appropriate 
for subsonic source motion. This is needed for ultimate 
implementation and numerical comparison of the dif- 
ferent formulations. 

A slightly modified integral formulation for the 
FW-H equation is needed because the current prac- 
tice is to assume that the FW-H integration surface 
corresponds to the body and is impenetrable. Equa- 
tion (6) is the appropriate form of the FW-H equation 
to start the development of an integral representation 
which has the same form as the traditional application 
of the FW-H equation. Following di Francesantonio, 14 
we define new variables Ui and Li as 

Ui = (1-A^ + — (14) 

Po Po 

and 

Li — PijTij + puiiufi r? n ) . (15) 

We have chosen a slightly different, but equivalent, 
definitions from that of reference 14 because p and pui 
are conservation variables often utilized in CFD codes. 
With these definitions, the FW-H equation may be 
written in its standard differential form: 

- A M(/)] + L [{PoUnW )] . (i6) 

This equation is particularly useful because Farassat’s 
formulation 1 A 6, 16 can be utilized directly to write an 
integral representation of the solution as 
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and j9g(x,t) can be determined by any method cur- 
rently available (e.g., see reference 8). In equation (17) 
the dot over a variable implies source time differentia- 
tion of that variable, Lm = LiMi , and a subscript r or 
n indicates a dot product of the vector with the unit 
vector in the radiation direction r or the unit vector 
in the surface normal direction n, respectively. 

Current rotor noise prediction codes can easily be 
modified to accommodate this new implementation of 
the FW-H equation. The major difference is that the 
integration surface is no longer restricted to the rotor 
blade surface and in addition to p', the values of p, 
pUi are needed as input. When the surface does cor- 
respond to the blade surface, the separation of source 
terms into thickness, loading, and quadrupole noise 
still has physical meaning; otherwise, the separation 
of the source terms into p' T , p ' L , and p'q is only math- 
ematical. Hence, the ability to give physical interpre- 
tation to the source terms continues to be a distinct 
and unique advantage of the FW-H equation. 


veloped to test the numerical implementation of equa- 
tion (17) without the quadrupole source term. The 
modified code is called FW-H/RKIR in this paper. 
RKIR was chosen as the platform to test the new 
FW-H implementation primarily because it already 
performs integration on a surface some distance from 
a rotor blade and has been coupled to the full poten- 
tial flow solver FPRBVI. 1819 A third code, WOP- 
WOP+, 8 which utilizes the traditional FW-H imple- 
mentation (surface integration on the blade surface 
and an approximate quadrupole implementation) will 
also be used in the comparison. 



Numerical Comparison of Formulations 

Although we have shown analytically that the 
FW-H formulation has advantages over the Kirchhoff 
formulation, what really matters is how they com- 
pare in practice. Some comparisons have already been 
made (e.g., see references 5,14 and 17). In reference 14, 
di Francesantonio concluded that the main advantage 
of the FW-H equation applied on a Kirchhoff-type in- 
tegration surface is that interaction with CFD codes 
is easier because the normal derivative of pressure is 
no longer required. If this is the only advantage, and 
indeed we recognize that the normal derivative calcu- 
lation can be cumbersome, a simple solution would be 
to make the substitution 


dp ^ dpUi 

d n = ~ Ui ~dT 


(18) 


in equation (8). This result is just the linear momen- 
tum equation, which is applicable in the linear flow 
region. Nevertheless, we believe there are other ad- 
vantages which we will now demonstrate numerically. 

For this work, a new computer code based on a mod- 
ification of the RKIR code (Rotating KIRchhoff for- 
mulation) developed by Lyrintzis et al. 10 has been de- 
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Figure 1. Comparison of predicted and measured 20 
acoustic pressure at an in-plane observer location, 
3AR from the rotor hub of an untwisted UH-1H model 
rotor in hover ( Mh = 0.88). 

The first comparison is for an untwisted UH-1H 
model-scale rotor operating in hover with a hover-tip 
Mach number Mh = 0.88. 20 Figure 1 shows a com- 
parison of acoustic pressure time history for both the 
Kirchhoff and FW-H methods on a integration surface 
which was located approximately 1.37 chords away 
from the rotor in the direction normal to the blade 
surface and extending 1.25 chords beyond the blade 
tip. The full potential computation was performed on 
a 80 x 36 x 24 grid, which is somewhat coarse. The 
two computations are almost indistinguishable in this 
case — an indication that the integration surface is in- 
deed in the linear flow region. The underprediction 
of the negative peak is a result of using a coarse grid. 
Brentner et al. 5 found that the agreement is improved 
with a finer grid. Small oscillations in the signal, near 
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Figure 2. Cross section showing the location of the 
integration surfaces with respect to the rotor blade. 
The vertical distance from the blade chord, in units of 
chord length, are labeled z/c. The value of the grid 
index normal to the blade is labeled k. 

the two positive peaks, are evident in both the Kirch- 
hoff and FW-H solutions. These oscillations are al- 
most certainly due to inaccurate quadrature over pan- 
els moving at high speed. The oscillations disappear 
as the integration surface size is reduced. 

Now that the FW-H/RKIR code has been intro- 
duced, we wish to examine the sensitivity of each for- 
mulation to the placement of the integration surface. 
Brentner et al. 5 found that the Kirchhoff solution var- 
ied somewhat with location of the integration surface. 
Figure 2 shows a cross section of five different inte- 
gration (Kirchhoff) surface locations ranging from one 
grid line off the surface to 1.37 chordlengths off the 
surface. The Kirchhoff acoustic pressure predictions 
from RKIR code for each of these surface locations 
are shown in figure 3. As the integration surface is 
brought nearer to the surface and the input data is 
no longer compatible with the linear wave propaga- 
tion assumption, the predicted acoustic pressure be- 
comes meaningless. Although expected, this aspect of 
the Kirchhoff method is troublesome. If the surface is 
not positioned properly the error can be substantial. 
Worse yet, if the integration surface is just positioned 
slightly in the nonlinear region the solution may be 
significantly in error but not enough so to be easily 
recognized. 

Figure 4 shows the noise prediction using the FW-H 


Figure 3. Comparison of predicted acoustic pressure 
using the Kirchhoff formulation with varying integra- 
tion surface locations. These predictions are for an 
observer located 3AR from a UH-1H model rotor hov- 
ering at Mh = 0.88. The experimental data is from 
reference 20. 



Figure 4. Comparison of predicted acoustic pressure 
using the FW-H formulation integration surface lo- 
cations. These predictions are for an observer lo- 
cated 3AR from a UH-1H model rotor hovering at 
Mh = 0.88. The experimental data is from reference 
20 . 
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Figure 5. Comparison of noise components predicted 
by the FW-H/RKIR and WOPWOP+ codes for a 
hover UH-1H model rotor ( Mh = 0.88, inplane ob- 
server 3AR from rotor hub). 

formulation given in equation (17) for the same set 
of integration surfaces and CFD input data as shown 
in figure 3. The volume quadrupole source, which 
exists only outside the integration surface, has been 
neglected in this calculation. The advantage of the 
FW-H formulation is clear: for an integration surface 
near or on the physical body, the predicted acoustic 
signal is essentially that of thickness and loading noise 
alone. As the integration surface is moved farther and 
farther away, more and more of the quadrupole source 
contribution is accounted for by the surface integrals. 
Hence, we would say that the principal advantage of 
the FW-H formulation for aeroacoustics is the relax- 
ation of integration surface placement restrictions. In 
fact when the volume quadrupole source is included in 
the noise computation, the location of the integration 
surface is only a matter of choice and convenience. 

Another traditional advantage of the FW-H method 
is the physical basis and identification of the source 
terms. If equation (17) is used on a surface away from 
the body this feature is not retained, however, a second 
computation can be made on the body surface to deter- 
mine thickness and loading noise. This has been done 
in figure 5, which is a comparison of FW-H/RKIR 
predictions with a WOPWOP+ prediction. Two FW- 
H/RKIR computations are show in figure 5: an inte- 
gration surface coincident with the rotor blade surface 
to predict thickness and loading noise, and an inte- 


Figure 6. Schematic showing three inplane microphone 
locations used in the the measurement of noise from 
the model scale Operational Loads Survey (OLS) ro- 
tor. 21 


gration surface located approximately 1.5 chordlengths 
away from the blade to predict the total noise. Note 
that the thickness noise predictions from WOP WOP + 
and FW-H/RKIR are identical and there is only a 
small difference in the predicted loading noise. The 
difference in the predicted loading noise is due to a 
difference in how the integration over the blade tip 
face is handled. The total noise, which includes the ef- 
fect of the quadrupole, is also in very close agreement 
even though the volume used in WOP WOP + is not 
identical to the region enclosed in the FW-H/RKIR 
surface integration. The negative peak is also in bet- 
ter agreement than the earlier figures because an Euler 
solution from Baeder 17 was used as input rather than 
the FPRBVI solution used in figures 1, 3, and 4. 

A model-scale test of the Operational Loads Sur- 
vey (OLS) rotor is selected for a final comparison. 
The predicted noise from FW-H/RKIR, RKIR, and 
WOPWOP+ are compared with experimental data 21 
at three inplane microphone positions, shown schemat- 
ically in figure 6. The rotor was operating in a for- 
ward flight condition with advancing-tip Mach number 
Mat — 0.84 and advance ratio /i = 0.27. A FPRBVI 
solution (80 x 36 x 24 grid) was used as input data 
for all three noise predictions shown in figure 7. All 
of these predictions agree quite well with the data — 
both in directivity and amplitude. All of the codes un- 
derpredict the negative peak pressure for microphone 
6, but this is most likely attributed to the FPRBVI 
solution rather that the noise prediction codes. The 
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Figure 7. Comparison of predicted and measured 21 acoustic pressure at three microphone locations for the model 
scale Operation Loads Survey (OLS) rotor ( Mat = 0.84, p — 0.27). 


differences between the predictions is most noticeable 
in the positive peaks, but even there predictions vary 
from each other by no more than 10 Pascals. 

A New Metric for Comparison 

The question of where to place the Kirchhoff surface, 
and the analogous question of how far out to perform 
the quadrupole integration, have a strong impact on 
the decision of which method is most efficient com- 
putationally. Two things must be considered: i) the 
accuracy and extent of the CFD computation needed 
as input data for noise prediction, and ii) the amount 
of input data required. Both of these will depend on 
the size of the nonlinear region surrounding the body 
generating the noise. While at present we are unable 
to give a completely satisfactory answer to the above 
questions, we can provide some guidelines using a new 
metric as follows. 

We note that the solution of the FW-H equation 
with the quadrupole source term invariably involves 
Tij and its first and second derivatives. Therefore, 
it is imperative that not only T \j is calculated accu- 
rately, but also its first and second derivatives in the 
source region. Similarly, the Kirchhoff formula tells us 
that on the Kirchhoff surface, p' , p 1 , and dp 1 /dn = p' n 
must be computed accurately in the CFD solution. 
This indicates that the error analysis in all high resolu- 
tion CFD computations must be based on the Sobolev 
norm. This norm is used very often in finite element 
analysis 22 and we propose such a norm in aeroacous- 
tics. We will not present any numerical results in this 
paper based on the Sobolev norm. 

We first address the problem of how to compare two 
high resolution CFD solutions for both FW-H and 


Kirchhoff methods. We assume all parameters of CFD 
and acoustic calculations are dimensionless. Let V be 
the volume where CFD computations are performed 
with the boundary dV. We define the Sobolev norm 
of j as 
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where T is a convenient time period usually taken as 
the inverse of the blade passage frequency. We have 
now a metric for comparing two CFD calculations as 
follows. We define the distance (error) between two 
results by 


d(Tij,T?j) = \\Tlj — ( 20 ) 

where T\- and T- 2 pertain to the two sets of compu- 
tational results. We may agree that the error is small 
if 

II'tU rp 2 II 

11 |j Tl || k « 1 • (21) 

The two sets of results may come from two different 
CFD computations. 

Now we consider the Kirchhoff method. Assume S 
is the Kirchhoff surface over which the nondimensional 
p', p', and p r n are specified. We define a Sobolev norm 
of p 1 and distance for two solutions from CFD compu- 
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tations p' 1 and p' 2 as follows: 


in 

0 S' 


\P '\ 2 + IP? + \P'U I 2 


'I 1/2 

dSdt\ 


( 22 ) 


d(p' 1 ,p' 2 ) = \\p n -p' 2 \\ s ■ (23) 

We can use this distance or error function to know 
when to stop a CFD grid refinement. Unfortunately, 
this norm would not tell us when we are in the linear 
region or whether the dispersion and dissipation errors 
have substantially influenced p', p', and p' n . These ef- 
fects are governed by grid size as well as artificial vis- 
cosity. In the study of these effects in high resolution 
CFD calculations, we must employ a Sobolev norm in 
defining the computational errors. 

An alternate use of the norm defined in equa- 
tion (19) is to decide the volume of quadrupole source 
included in our noise calculations. Let V\ and V 2 be 
two volumes such that V\ C U- Then assume that 
T\- — 0 outside Vi. Using the Sobolev norm with vol- 
ume integration over V 2 , we can say that Vi includes 
all quadrupoles needed for noise calculation if 


rpl rp 2 

± ij ± ij V 2 

T T. 

Vi 


< 1 . 


This means that 


*jllV 2 \Vi 




T 1 . 

Vl 


(24) 


(25) 


where V 2 \ Vi is the volume enclosed between dV\ and 
dV^- This answers how far from the blade surface we 
must include quadrupole sources. 


Conclusions 

In this paper we have compared two useful aeroa- 
coustic tools: i) the Lighthill acoustic analogy as em- 
bodied in the FW-H equation, and ii) the Kirch- 
hoff formulation for moving surfaces. Both of these 
methodologies have proven their usefulness in rotor 
noise prediction. It is because both methods work well 
that it is difficult to determine which to use for a par- 
ticular application. In a comparison of the governing 
equations, we have shown that the FW-H approach 
can include nonlinear flow effects in the surface in- 
tegration if the usual assumption of an impenetrable 
surface is relaxed. In fact, when the integration sur- 
face for each is in the linear flow region we have shown 
that the FW-H equation is equivalent to the Kirchhoff 
governing equation. 

The FW-H equation is based on the conservation 
laws of fluid mechanics rather than the wave equation, 
as is the case for the Kirchhoff formula. As a conse- 
quence, the FW-H equation is not appropriate for all 


types of wave propagation (e.g., the FW-H equation 
is not appropriate for electromagnetic wave propaga- 
tion, while the Kirchhoff formula could be utilized). 
But the superiority of the FW-H for the aeroacoustics 
of rotating blades has been demonstrated through sev- 
eral numerical examples in this paper. The placement 
of the integration surface is a matter of convenience as 
long as the quadrupole source is utilized. The FW-H 
method also has the advantage that it separates the 
predicted noise into physical components (i.e., thick- 
ness, loading, and quadrupole), explicitly. The Kirch- 
hoff method does not offer this insight into the nature 
of the acoustic field. 

It is well known that the quadrupole sources are 
responsible for noise generation as well as distortion 
of the acoustic waveform. The intense quadrupole 
sources are in the vicinity of the blades. Therefore, 
if we use a surface which encloses the blade and the 
volume of intense quadrupoles in the FW-H method, 
we can calculate the level of the acoustic pressure ac- 
curately. The role of the weaker quadrupoles, which 
are farther away from the physical body, is primarily 
to provide a small distortion to the acoustic waveform. 
Hence, even when the integration surface is fairly close 
to the noise generating surface, it may be acceptable 
to neglect the external quadrupoles. In comparison, 
the Kirchhoff formula can predict acoustic pressures 
that are substantially in error if the Kirchhoff surface 
is located inside the nonlinear region. The nature and 
order of magnitude of this error may be hard to esti- 
mate or even recognize. 
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